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DETERMINATION OF THE RATIO OF ABSORPTION COEFFICIENTS AT DIFFERENT WAVELENGTHS IN A 
SCATTERING MEDIUM 

This invention relates to the determination of the ratio of absorption coefficients at 
different wavelengths in a scattering medium, and to the determination of the relative 
concentrations of chromophores in the medium. 

There are a number of industrial processes in which measurement of the ratio of 
chrornophore concentrations in a scattering medium would (if available) be an important 
indicator for process quality control. Examples are the measurement of water, protein and 
fat content in liquids and solids in bioreactors and food processing, and the measurement of 
relative dye concentrations in dyeing processes, which relative concentrations are critical to 
the consistency of the final colour achieved in the dyed product. 

Other industrial applications would be the measurement of the ratio of chrornophore 
concentrations in tissue samples, for example relative oxy-and deoxy haemoglobin 
concentrations. Measurement of such quantities in vivo would also be useful in support of 
clinical practice, as would measurement of the oxygen saturation of blood in tissue. 

It is an object of at least the preferred embodiment of the invention to provide a 
method of determining relative concentrations of chromophores in a scattering medium, and 
also to provide a method of determining relative absorption coefficients of the medium itself, 
since this may be sufficient for some applications. 

In a first aspect the invention provides a method of determining the ratio of the 
absorption coefficients and ^(X^ of a light scattering medium at two different 

wavelengths comprising 

a) passing light (as herein defined) through the medium 

b) measuring a characteristic X of the light at two wavelengths Xj and Xj w ^ch is 
affected by the absorption coefficient during its passage through the medium 

c) effecting changes Ajji^Xj) and Ajx^Xj) in the absorption coefficient at the two 
wavelengths; 

d) repeating step (b) 

e) determining the ratio of the absorption coefficients lij^x) anc * M-^XJ at said two 
wavelengths from said changes therein or from quantities representative thereof and from the 
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characteristic X and from the transport scattering coefficients of the medium fi^i) and fi/^ 
at the two wavelengths. 

The term "light" as used herein includes the near infrared spectrum. 

The said changes in the absorption coefficient may be effected by a common cause, 
for example by varying the concentration of a chromophore in the scattering medium, the 
quantities representative of the changes in the absorption coefficient being the absorption 
coefficients Afx^), A^(^z)- In tissue, these changes can be deliberately generated by 
injecting an exogenous dye (e.g. Indocyanine green) intravenously or modifying haemoglobin 
oxygenation or volume. 

The method may comprise first adding to the medium the chromophore where 
concentration is varied. 

The ratio of the absorption coefficients a "d H a (X2) may be determined from the 

ratio of the chromophore absorption coefficients Afx a (^i)> Ap^X^) and from the ratio of the 
transport scattering coefficients of the scattering medium ^'(XJ and ^Q^)- 

In one embodiment the characteristic measured in step (b) may be the attenuation A 
of the light during its passage through the medium. 

In another embodiment the light may be intensity modulated or pulsed light, the 
measured characteristic X being the phase or depth of the modulation or the time of flight 
(transmission time) of the light through the medium respectively. 

The invention also provides a method of determining the variation with wavelength 
of the absorption coefficient of a scattering medium (the spectral characteristic) comprising 
performing the method as set forth above at a number of wavelengths across a waveband of 
interest, and in step (e) thereof determining the ratio of the absorption coefficients relative to 
that at a reference wavelength the variation of said ratios with wavelength being the said 
spectral characteristic. 

In a further aspect the invention provides a method of determining the ratio of 
chromophore concentrations of known spectral characteristics (as herein defined) in a light 
scattering medium comprising performing the method of steps (a) to (d) as set forth above 
at not less than n different wavelengths where n is the number of chromophores, one of said 
wavelengths being the reference wavelength X, 9 performing step (e) to determine the ratio 
M-aC^iV^aC^x) of the medium for each of the other wavelengths, and determining the relative 
concentrations of the chromophores in the medium from said ratios. 
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By spectral characteristic we mean the variation of \i % with X. 

Where the absolute value of the absorption coefficient of one of the chromophores is 
known (eg water) or can be found by other means, the method will yield the absolute spectral 
characteristics and hence absolute concentrations of the other chromophores. 

The invention also includes apparatus for use in the method as set forth above. 

The invention will now be described merely by way of example with reference to the 
accompanying figures, wherein 

Figure 1 shows the fractional error in the ratio M^C^-O/^C^) from the approximation 
in equations 3 and 4, discussed below. Specifically, it shows fractional error in ratio 
^.(^/^aC^i) from approximations of equations 3 and 4 assuming \J^(Ki) = 0-02 mm" 1 , V^Q^) 
= fi s '(Xi) = 1 mm" 1 and a source detector r = 50 mm. 

Figure 2a shows the optical properties of a liquid tissue-simulating phantom used in 
an experiment described hereafter. The absorption coefficient is that of water (\i KW ) and the 
added dye (^a). The transport scattering coefficient (fx, 1 ) is based upon calculations from Mie 
theory for the concentration c, = 2% and the diameter (1 \xm) of the spheres. 

Figure 2b shows the spectrum of the absorption coefficient of the phantom normalised 
at X.=800 nm. Specifically, it shows spectrum of the relative absorption coefficient 
(^(X)/|ii 4 (800nm)) of the phantom normalised at X = 800nm. The dotted line shows the true 
values derived from the water and dye absorption of the phantom (n^ w + fi M in figure 2a). 
The measured attenuation changes that were induced by a small increase in \i Kd were used to 
estimate ^(X)/ji a (800nm) by assuming a constant \ij (dashed line, calculated from Eq. 6 
below) and by including the wavelength dependence of (V (solid line, Eq. 7). 

Figure 3 shows values of the transport scattering coefficient ^ measured on the heads 
of seven volunteers. The values are an estimate from measured TPSFs (data from Matcher 
et ai 1996). 

Figure 4 shows values of \ij measured on the calf muscle of eleven volunteers; |V 
measured on the head of 7 volunteers. The values are an estimate from measured TPSFs 
(data from Matcher et al 1996). 
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Figure 5 shows apparatus according to the invention. 

Changes in the chromophore concentration ie. changes in the absorption coefficient 
(A^iJ in a light scattering medium alter the intensity of the diffusely transmitted light 
collected at the surface of the medium. Here we describe how the spectral characteristic of 
the tissue absorption coefficient can be calculated from the known spectral characteristic of 
an injected dye absorption coefficient and the change in the measured spectrum of the 
intensity (attenuation), phase shift or mean time. 

Once the spectral characteristic of the tissue absorption coefficient has been 
determined, the relative concentrations of its chromophores, notably the concentrations of 
oxy- and deoxy - haemoglobin (Hb0 2 and Hb), can be calculated. An application is the 
monitoring of the oxygen saturation (Hb0 2 /(Hb0 2 + Hb)) of blood in tissue. 



Diffusion Theory Model 

For a semi-infinite halfspace the reflectance at a distance r from the light source can 



be written as 



P 2k- p 1 (1) 



In Eq. 1 p = (r l +zJ 2 y A and z 0 = l/u s ', where p s ' is the transport scattering coefficient. The 
velocity of light in the medium c = c,/n (where the speed of light in vacuum is c 0 ). 

<"*ff = (3-/*, C"»+y"»')) v4 is the effective attenuation coefficient and D = (S.^+jz,'))' 1 is the 
diffusion coefficient. The derivative of the attenuation A (defined as the logarithm of the 
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ratio of incident and detected intensity) with respect to change in is: 



dA 



d[L t 2-lnl0 l/p+M-ctr 



£ (2H. +»!.•). 



(2) 



Here it is suggested thai measurements of vAJaufJ^ can be used to estimate the product of 
absorption and scattering coefficient. For a known wavelength dependence of the spectral 
shape of can then be calculated from measurements of dAJfi x at different wavelengths. 



Assuming that scattering dominates absorption (ju s '» Eq. 2 can be approximated to: 

l2 



3p-iy . 

2-lnlO 



dA j _l 



1 

3^ s " 



(3) 



Furthermore, for large source-detector distances it is: 



(4) 



2-lnlO 



'_9A> 



-1 



■3^ $ * 



Eq. 4 states that the ratio of the absorption coefficients at different wavelengths (Xj, Xj) can 
be estimated from measurements of. attenuation changes (AA) at these wavelengths, the ratio 
of the changes in absorption coefficient and the wavelength dependence of /u^. 



AA(X 2 ) Ap-.fln) 
AA^)" A|l a (X 2 ) 



(5) 



To estimate the validity of the approximations leading from Eq. 2 to Eq.s 3 and 4, the errors 
in the calculated ratio MtQ^^Qi) are plotted in Fig. 1 as a function of the true Pi(kj)/fi&d 
assuming u^Xj) = 0.002 mm" 1 , fi,'Q*z) = Mm'Q^O = 1 mm" 1 and a source-detector distance r = 
50 mm. In the figure the fractional errors of Eq.s 3 and 4 are shown in solid and dotted line, 
respectively. The fractional errors increase with increasing differences of /^(Xj) and /^(X,). 
Variations in the magnitude of fx x of a factor of 2 result in errors of approximately 10%. 
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Determination of the relative absorption spectrum of a light-scattering phantom 

The method suggested above was tested in a liquid, light scattering tissue simulating phantom 
of known optical properties. Silica spheres (Monospher M1000, Merck, Germany) of 1 fan 
diameter were suspended in water and served as scattering centres in the phantom. The 
scattering properties could be derived from Mie theory calculations as the diameter of the 
spheres and the refractive index of both the spheres and the surrounding medium are known. 
The calculated scattering cross section of the spheres was compared with measurements from 
a collimated beam setup. Experimental and theoretical scattering cross sections were in 
agreement when the wavelength dependence of the refractive index of fused silica (n, = 1.447 
+ 3890X" 2 - 37X" 4 with X in nra) and water (n, = 1.32 + 6878X" 2 - 1.132.10" 9 X 4 ) were 
employed. The transport scattering coefficient £ was calculated for wavelengths between 600 
and 1050 nm for these refractive indices and the concentration of the spheres in the phantom 
(c, = 2% voL/vol.). The calculated decreases approximately linearly with wavelength and 
a first order regression through the calculated values gave /V=0.92 - 0.00047 Jifmm' 1 ] (X in 
nm). The absorption coefficient of the phantom consisted of two components: the water 
absorption (u^) and the absorption of a dye (S 109564, ICI, GB) (pj) added at a 
concentration of c d = 0.0031% volVvol.. The spectra of both^,' and/:, are shown in Fig. 2a. 



A CCD spectrometer in connection with a halogen light source was employed to measure the 
reflected light intensity of the phantom (volume 300 ml). The ends of the optical fibres that 
transported the light from the lamp to the phantom and from the phantom to the detector were 
submerged approximately 2 mm in the phantom at a source-detector distance of r = 38 mm. 
Two spectra (I 0 and I) of the reflected light intensity were recorded before and after the 
concentration of the dye was changed by Ac„ = 2.5 x 10 _4 %. This corresponds to a change 
in fi x of between 4% and 8% depending of wavelength. The spectrum of the attenuation 
change AA(X)=log[I 0 (X)/I(X)] was calculated. Ignoring the wavelength dependence of 
scattering gives a first approximation of the spectral characteristic /^(X) according to Eq. 4/5, 
normalised a) for the relative magnitude of the absorption change (ie. divided by the 
absorption spectrum of the dye pj) and b) for the value at X = 800 nm: 



MM 



AA(8Q0nm) H,. d (A.) 



MSOOnm) * l&M*-) Ji a . d (800nm) 



2 



(6) 
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This first approximation of /^(X)(dashed line in Fig. 2b) shows considerable discrepancies 
with the true ^(Xy^SOOnrn) spectrum (dotted line). Correcting the experimental data for the 
wavelength dependent scattering coefficient of the phantom (u,'(X) = 0.92 mm" 1 - 
0.00047Xfmm- 1 l. 



}jL a (800nm) 



AA(800nm) M- a .d^-) 
AA(X) V a ,d(800nm) 



^i s ' (800nm) 



(7) 



gives the spectrum shown plotted with a solid line in Fig. 2b. The experimental data are in 
good agreement with the expected fx x spectrum for °k < 950 mm. For longer wavelengths the 
experimental data are noisy and underestimate fx^ The poor signal to noise ratio for these 
wavelengths is due to the low sensitivity of the CCD detector and the higher absorption 
coefficient of the phantom. The likely explanation for the underreading at these wavelengths 
is, as demonstrated in Fig.l, that the approximations used for the derivation of Eq. 3-5 break 
down for larger ratios ^(^//^(SOOnm). 



Calculation of relative cbromophore concentrations from the relative absorption 
spectrum 

In the experiment, the determination of ratio of chromophore concentrations and the 
oxygenisaturation of blood depends on the following points: 

a) Attenuation spectra can be corrected for the wavelength dependence of fjQ\ Figure 3 
and 4 show fi* of the head and the calf muscle from volunteers obtained from 
measured temporal point spread functions (TPSFs). The absolute ju s * varies by up to 
30% for different volunteers, however it has a similar wavelength^ characteristic. 
Therefore, a correction of the attenuation spectra for variations is feasible. 

b) A pure change in absorption coefficient occurs and the wavelength dependence of the 
absorption coefficient change is known. Absorption changes can be induced by an 
exogenous dye or by changes of haemoglobin concentration. Change in Hb or Hb0 2 
can be generated by small changes in the inspired gas. Alternatively, the absorption 
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coefficient can be reduced by dilution (for example by injecting saline). 

c) The spectral characteristic of all tissue chromophores is known: 

For tissue and wavelengths in the near infrared the dominant chromophores are 
deoxy- and oxyhemoglobin and water, the absorption spectra of which are known 
In general the absorption coefficient at a wavelength X can be written as 



where s a (unit: mm" W) is the extinction coefficient of the n* chromophore and 
c a is ,ts concentration (unit: mM). Once the ratio ^)/*(\) at least at, n 
wavelengths has been determined ( K is a reference wavelength), an iterative technique 
can be used to calculate the relative chromophore concentrations c^ (where c, 
signifies the concentration of a reference chromophore). 

The number of wavelengths has to be equal to or higher than the number of chromophores 
m the medium. The absorption spectrum of the dye used to induce attenuation changes has 
to be known. For the determination of the relative oxy- /deoxy haemoglobin concentrations 
in tissue, the dye could for example be Indocyanine green. 

In the discussion above measured changes in the light intensity were used. Equally changes 
in the time of flight, the phase shift or the modulation depth of an intensity modulated light 
wave can be exploited. 

Thus changes in the mean (transit) time (<t>) of the photons or the phase (*) of an intensity 
lightwave (modulated at the frequency Vm ) can be used to compute the spectral characteristic 
of n,. 

The mean transit time is given by 

<t> = BL (8) 

2c-(D+p-7tvD) 

and the phase is approximately a linear function of <t>: 
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cD = -2tc-v m - < t > 



The derivation of <t> with respect to fx, can be approximated as 



8<t> = ~ p JT. £Z._L 
3jj. a 4c % ~ ^fi a }i a 



Therefore, 



a <t> ^ 



-2 



lV3 



(ID 



The equivalent to Eq.5 with respect to measured changes in mean transit time is therefore 



A<t>(Xo) A|X a (X t ) 



_A<c>(X l ) A(X a (X 2 ) 



(12) 



The measurements of phase changes AO, Eq. 12 can be combined with Eq. 9 to give 



JA*^) Afi a U 2 ) 



(13) 



The modulation depth can be written as 



where ip 0 = n 6ff .p[l + X 2 ] ,/4 , % = "^o -sin (6/2), 
■4), = o|; 0 .cos (6/2), 6 = tan -1 (x). •*!><» = IWP 31111 
X = (23tVM)/(fA 4 c). 

Changes in fi, induce changes AM. It can be shown that equivalently to Eq.7, 12 and 13 the 
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ratio of absorption coefficients at two different wavelengths can be derived from the measured 
fractional change 6M = AM/M at these wavelengths according to 



6MU 2 ) .An.UJ 



[[6M(A X ) Au a U 2 ) 



H*U 2 ) j 



l/s 



Figure 5 shows diagrammatically apparatus according to one embodiment of the 
invention in which changes in attenuation of light passing through the medium are measured. 
A medium 10 under examination receives light from a white light source 12 via an optical 
fibre 14. 

Diffused light is received by an optical fibre, and is taken to a wavelength dispersive 
device 18 and thence to a detector 20. The output of the detector 20 is fed to a PC or other 
data processor 22. 

The wavelength dispersive device, for example a spectrometer comprising an optical 
slit and a diffraction grating divides up the white light issuing from the medium into its 
constituent wavelengths, those of interest then being selected by output slits or by being 
detected by respective elements if the detector 20 is of the staring-array type. The detector 
20 may be a charge-coupled device or an array of photodiodes. 

Provided the intensity of the light source 12 is kept constant, it is unnecessary to 
provide the detector 20 with a measure of the input light intensity. All necessary data is 
contained in changes in intensity arising from imposed changes in X and iv 

As an alternative to the spectrometer, a wavelength selective device may be employed 
such as a rotatably mounted filter wheel with interference filters with two or more coloured 
filter elements for passing the wavelengths of interest. 

The wavelengths of interest are typically in the range 650nm to lOOOnm'for tests on 
tissue. When the distance between light source and detector fibers is small (<5mm) 
wavelengths of up to about 2000nm (2u.m) are appropriate. Preferably the wavelength 
dispersive or selective device should have a bandwidth of less than 10 nm. 

In an alternative embodiment the white light source 12 is replaced by switchable or 
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otherwise selectable monochromatic light sources (eg laser diodes or LEDs) of different 
wavelengths, the wavelength dispersive or selective device 18 then not being required. Laser 
diodes are preferred to LEDS because the wavelength of the emitted light is more 
monochromatic. 

TIM *U« tr\ Vio l'mnlpmPnfpH hu mP5Cimn(7 thp. time nf fHpht <t>. a DUlSed 

VYliC>il Llllx UlVbUUUIl AO IW iAliy»wiiiv*4kvv* */; *** — » — — cj * * 

laser diode is used as the source 12 to provide pulsed light of a specific wavelength to the 
medium 10. 

The wavelength dispersive device 18 is omitted and the detector 20 is replaced by 
means to determine the time of arrival of the pulse. A signal is provided via line 24 from 
the drive circuitry of the laser diode to indicate the time of initiation of the pulse to the 
detector so that <t> may be established. 

For implementing the invention by measuring the phase shift A<|> or depth of 
modulation of Light, light source 12 is an intensity modulated monochromatic source with a 
choice of frequencies, either switchable laser diodes or LEDs, or a white light source with 
several alternative filters. The detector 20 is phase-sensitive, and receives a reference signal 
from the input to the sample via a line 26. Phase detection may conveniently be by a lock-in 
technique in which the output from the medium is compared with the input from line 26 by 
accurately multiplying the detector signal by both the reference signal and the reference signal 
shifted in phase by 90°. 

Each feature disclosed in this specification (which term includes the claims) and/or 
shown in the drawings may be incorporated in the invention independently of other disclosed 
and/or illustrated features. 
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CLAIMS 



1. A method of determining the ratio of the absolute absorption coefficients and u,^) 
of a light scattering medium at two different wavelengths comprising 

a) passing light (as herein defined) through the medium 

b) measuring a characteristic X of the light at two wavelengths X, and X, which is affected 
by the absorption coefficient during its passage through the medium 

c) effecting changes Aa,(^) and Au,(^ in the absorption coefficients at the two 
wavelengths; 

d) repeating step (b) 

e) determining the ratio of the absolute absorption coefficients u,^ and at said two 
wavelengths from the said changes therein Af^X,) and Af^ or from quantities 
representative thereof, from the characteristic X, and from the ratio of transport scattering 
coefficients of the medium fV(X,) and u.^ at the two wavelengths. 

2. A method as claimed in claim 1 comprising effecting said changes in the absorption 
coefficients by means of a common cause. 

3. A method as claimed in claim 1 or claim 2 comprising effecting said changes in the 
absorption coefficient A^XJ and Au^ by varying the concentration of an existing 
chromophore in the scattering medium. 

4. A method as claimed in claim 1 or claim 2 comprising effecting said changes in the 
absorption coefficients Au^X,) and Au^ by adding a chromophore" to the scattering 
medium. 

5. A method as claimed in claim 3 or claim 4 wherein the ratio of the effected absorption 
coefficient changes Au^X,) and Au^ is determined from the ratio of the absolute 
chromophore absorption coefficients fx ld (Xi), \i lA Q^). 
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6. A method as claimed in any preceding claim wherein the characteristic X measured in step 
(b) is the attenuation A of the light during its passage through the medium. 

7. A method as claimed in any of claims 1 to 5 wherein the light is intensity modulated light, 

tli* moocnrp^ nUomntf>Ti<:tir> V Vainer thr nhacp ar intp-ncitv nr Hpnth nf thft modulation nf the 

light through the medium. 

8. A method as claimed in claims 1 to 5 wherein the light is a short pulse of light, the 
measure characteristic X being the time of flight (as mean time of transmission) of the light 
through the medium. 

9. A method of determining the variation with wavelength of the ratio of absorption 
coefficients of a scattering medium (the relative spectral characteristic) comprising performing 
the method of any of claims 1 to 8 at a number of wavelengths across a waveband of interest, 
and in step (e) thereof determining the ratio of the absorption coefficients relative to that at 
a reference wavelength the variation of said ratios with wavelength being the said spectral 
characteristic. 

10. A method of determining the ratio of concentrations of chromophores of known spectral 
characteristic (as herein defined) in a light scattering medium comprising performing the 
method of steps (a) to (d) of any preceding claim at not less than n different wavelengths 
where n is the number of chromophores, one of said n wavelengths being the reference 
wavelength performing step (e) to determine the ratio \i z Q^i)/^(K) of tlie medium for each 
of the other wavelengths, and determining the ratio of concentrations of the chromophores in 
the medium from said ratios. 

11. A method of determining the absolute concentration of chromophores of known 
spectral characteristic, as claimed in claim 9, where one of the chromophores is of known 
concentration. 

12. Apparatus for determining the ratio of the absolute absorption coefficient \iJQ^ and 
P-a(^i) °* a ^gkt scattering medium at two different wavelengths comprising 
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a) means for passing light (as herein defined) through the medium; 

b) means for measuring a characteristic X of the light at least two wavelengths \ and X, 
which is affected by the absorption coefficient during its passage through the medium; 

c) means for effecting changes A^) and Au^ in the absorption coefficient at the two 
wavelengths; 

d) means for determining the ratio of the absolute absorption coefficients u^) and u,^ 
at said two wavelengths from the said changes therein Au,^) and Au,^ or from quantities 
representative thereof, from the characteristic X, and from the ratio of transport scattering 
coefficients of the medium ^(XJ and fx,'(KJ at the two wavelengths. 

13. Apparatus as claimed in claim 12 wherein the means (b) is adapted to measure the 
attenuation A of the light during its passage through the medium as said characteristic X. 

14. Apparatus as claimed in claim 12 wherein the light is intensity modulated light and 
means (b) is adapted to measure the phase or the intensity or the depth of modulation of the 
modulated light through the medium as said characteristic X. 

15. Apparatus as claimed in claim 12 wherein the light is a short pulse of light, and means 
(b) is adapted to measure the time of flight (mean time of transmission) of the light through 
the medium as said characteristic X. 

16. Apparatus for determining the ratio of concentrations of chromophores of known spectral 
characteristic (as herein defined) in a light scattering medium comprising means performing 
the method of steps (a) to (d) of any of claims 1 to 11 at not less than n different wavelengths 
where n is the number of chromophores, one of said n wavelengths being the reference 
wavelength X„ means for performing step (e) to determine the ratio Au,^) and Aji.^ of 
the medium for each of the other wavelengths, and means for determining the ratio of 
concentrations of the chromophores in the medium from said ratios. 

17. Apparatus as claimed in any of claims 12 to 16 wherein the means (c) is arranged to 
effect said changes by a common cause. 
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